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ABSTRACT 


This study investigates fin-line filter structures, a 
type of E-plane waveguide device. Expressions germane to 
the analysis of fin-line structures are developed. A 
computer aided design and analysis program based upon a 
mode-matching technique is described. Filters designed 
by this program are fabricated and tested in X-band. Good 
agreement with predicted response is obtained. Filters 
fabricated and tested in Ku and Ka bands by other researchers 
are analyzed by the new program. Good agreement between 


predicted response and published performance is noted. 
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.I. INTRODUCTION 


Integrated circuits for use below about 3 GHz demonstrate 
clear advantages in terms of size, weight, cost and electri- 
cal performance, over alternative techniques. These 
advantages often accrue as a result of component miniaturi- 
zation. However, as the frequency is increased to the 
centimeter and millimeter regimes, the resultant excessive 
miniaturization creates problems associated with critical 
mechanical tolerances, and questionable production uniformity. 
Moreover, concomitant undesireable electrical performance 
involving radiation loss, spurious coupling, dispersion and 
higher order mode propagation becomes unacceptable [Ref. lj. 
Over the last decade, integrated fin-line structures have 
received increased attention as attractive media for 
low-insertion loss design applications which avoid excessive 
miniaturization, yet offer the potential for low-cost 
batch-processed production techniques. 

The fin-line filter structure consists simply of 
eemductive strips, or fins, which bridge the broad walls of 
a rectangular waveguide. That is, for waveguide operated 
ina TE o mode, the conducting fins are suspended in the 
E-plane. Metal sheet structures suspended in the center of 
the waveguide (herein referred to as "Class I" structures) 


offer the simplest possible structure for analysis and 


“24 
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ааа = 
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design. These unilateral center-line fins present fart 


mean and suspension difficulties. Bilateral Tin-lane Filters 
may be constructed as printed circuits on low-loss dielectric 
substrates and inserted in the center of the guide. These 
ШҰГавв ІЛ") filters offer the advantage of Construction by 
well developed pp cipeUwrt techniques , Vand ЕтуеТу 
straight forward analysis [Refs. 2, 3], but still present 

the problem of mechanically fixing the fin-line structure in 
the center of the guide. Clearly, fin-line structures printed 
on dielectric substrate and affixed to the waveguide side- 
walls ("Class III" structures) offer mechanical advantages 
Omer the other two class types. The authors know of no 
previous investigations concerning Class ITI strucutres. 

The purpose of this thesis study was to develop and 
validate a computer aided design (CAD) program for the design 
and analysis of all three class types cf E-plane fin-iine 
Filters. To this end, the appropriate mathematical 
expressions were derived for a generalized fin-line structure. 


The generalized structure has two different dielectric 
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| i} 
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Ttenacceenı ned above 


ct 


materials in the guide. The Class IT 
corresponds tc the case where the higher dielectric constant 


` = 


m Trial is in the center of the guide. The Class III 


Ih 


|-- 


ult 


(D 


E. 
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pa 
к! 


has the higher permittivity material along the waveguide wa 
The generalized Class I filter has the same permittivit; 


throughout the fin-line structure. 





The derived transcendental equations relating the mode 
wave number values in different zones of the waveguide 
apply to all class types. However, when the difference 
in dielectric constants is large the trigonometric functions 
in the transcendental relationship may become hyperbolic. 
Moreover, Class II and Class III structures require different 
but equivalent expressions in the numerical analysis to 
avoid singularity occurances. 

The analysis method employed is a mode-matching technique 
uniquely applied to include the effects due to an arbitrary 
number of higher order modes. A convergence dependence 
analysis assures the most efficient regional distribution of 
the limited number of higher order modes retained for the 
numerical analysis. 

Validation of the complete computer program for Class , 
шта. 11 filters was initially conducted by analyzing k, and 
к. band fin-line filters designed and constructed by Shih 
E. 2], and Arndt, et al LRef. 3].  FQnally,.several Class 
I and Class III filters were designed, constructed and 
tested locally. These test filters were resonant in 


X-band (8-12 GHz) to match available test hardware. 





TEARS 


For the types of fin-line structures in waveguides 
considered here, each geometrical discontinuity may be 


considered in two distinct coordinate frames. It is useful, 


e 


paucPerore, to expand the unknown fields in Beth regions in 
terms of their associated normal modes. Field continuity 
requirements are then imposed on the modes at the interface 
e" the regions. Finally, the orthogonality property cf che 
normal modes is used to derive an infinite set of linear 
Simultaneous equations for the unknown modal coefficients. 
This is the essence of the mode-matching technique for the 
solution of boundary value problems [Ref. à]. Moreover, 
this procedure leads directly to the scattering matrix for 
ie discontinuity. Since the fin-line is a series of such 


discontinuities separated by transmission line sections, a 


cascading technique may be employed to develop an equivalent 


rh 


seattering matrix for the entire structure. 


1-4 
rt 


fects of all the possible modes upon 


(0 
(0 
rh 


pa mel ple, tla 
each of the others is inherently included, and the analysis 

is exact. In practice, the number of modes considered must 

Bel limited to a reasonably finite set. Accurate and 

efficient convergence of the subsequent numerical calculations 


is achieved by following tne guidelines suggested by the 


КЕЕ св ог Shih and Gray [Rer. 51. 
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ене demiyatıom of the eauiwelent se@fterine matrix for 


Mh 


led 
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the ceneralized fin-line structure follows. Deta 


0Q 


intermediate steps of the development are contained in the 
appropriate appendices. The basic geometry of the E-plane 
Aline is shown in Fic. IL. It is important to note 

that the assumed "magnetic wall" (x=0), along which the 


parallel magnetic field is zero, implies TE G (m odd) 


Si 


"ТЕЛ 


ineident waveforms. The computer simulation is based upon a 


ТЕ) incident weve. 


ERO EXERESSIONS FOR NORMAL MODES AND EIGBNVALUES 


HF) 


By assuming incident modal electric (E) fields of the 


ES type, field symmetries allow the simplified waveguide 


= . 
zn 
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geometry shown in Fig. II-2. The conducting fin is assumed 


UG 


ل 


de 
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. The dđdielectrie materials im ue gu 
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are assumed to be lossless, homogeneous, and istropic. Also, 
the fin-iine structures and substrates must be constructed 
to maintain symmetry about the x = 0 plane. For simplicity 
of discussion, the guide is considered to be comprised of 

two regions, and two zones. The regions are Z dependent. 


uide before the metal fin boundary, 


Gal 
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el: 
ЕУ 
(D 


Region I refers 
and Region II refers to the guide after (plus Z direction) 


SNS 


el: 


the fin boundary. The zones are X dependen 
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sometimes calied the lower zone, refers to the condition: 
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while Zone 2, sometimes called the upper zone, refers to 


The condition: 


Modal eigenvalues in Region I use lower-case subsceripts, 
while those of Region II use upper-case. Regional divisions 
Ec occassionally indicated explicitly bv сае superseripes, 
шүн and де 


The 4 dependence of all fields is of the form: 


: х E ei 
where Y. 13 the propagation constant of кпе m mode. 


Likewise, the time dependence of all fields is of the form: 


pa 
ba 
ct 


where » 1s the radian frequency of the electromagnetic (EM) 
wave. These terms are understood to multiply all field 
expressions, and will not be explicitly shown, except where 
ÁS more convenient, or constructive to do lso 

The general solution to the wave equation has the form: 


ae t=. t 
[ae сое кк Ро о ИК Ј 
| J K JS J & 


4 


LE 





where the E are the associated wave numbers and the x, аге 
the appropriate directional variables, x or y. 

For our purposes, the electric field in the Ү есен 
E is of greatest interest. In Region I, the appropriate 


expressions for the normal modes of the 2 field are found 


to be (Appendix A): 
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The wave numbers for the u mode are related Dy the 
y 


transcendental equations: 





m C c um D А. 2 
in n O u % 

{ 

] Az ] = Ë a- ТА, 

Kin an К , 9 Kon CO Kn d) me) 

where 
k I 2T = 271 И 
O A С 
O 


Equation 11.4.2 is the usual dispersion relation, 
wnere с. апа e, are the redative dielectric CoOmStanes Om ene 
me zones. The transcendental equations (11.A.2) amd (II.À.5) 
are solved simultaneously to determine the value of the wave 
number for each mode in both zones. 


The corresponding expressions for Region II are similarly 


Oil tO be: 


(2317? COSI x ТІ Ф 
Y cu r d КЕ = == 
E s (ТЕЖА U) 
D асык 
CE суас _ (2-2) ea ya 
a-d UP = = 





= 
The wave number for the p | mode is: 


_ (2P-1)r : 
Кз, = > DOS 
E Tr . 
Kup * Та-аў ша. 


It should be noted that in Region II, the eigenvalues for 


Kip Fe uncoupled from the eigenvalues for Kup 


mode number, p, runs from 1 to e individually in each zone. 


so that the 


h . th І | 
Ihe propagation number for the m mode in Region I is 


given by: 


71 = OS СЕ e, Kk и 
m im 1 o 


The Region II propagation numbers are given by: 


Positive real values for the propagation number corresponds 
to evanescent waves. Propagating wave modes occur when 


the propagation number is imaginary. 


HE 





B. MODAL ANALYSIS 
Any physically realizable electric field may be written 
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ELEM оле lines combination or the normal modas., That is: 


ЕР = z A $. 
У me il 
and 
E ЕХ ӨШ : 
y Se D 


The y and ШЕ are the normal modes or the m meld in Region 
I and Region II, respectively. The i and D. are the 
associated amplitudes. 

The interface discontinuity may be analyzed as a two 
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Bere junction, with values as defined in 


boundary conditions on the $ апа Y, [Ref. 4] require that: 
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А. - Coefficients of Modes Ineident from the Left 

E ="Coctticients of Modes Jeaving Junction to the Left 
E cNCocséuejents of Meses Incident zen the Risch 

D. = Coeffiezents of Modes leaving junetson to the віват. 
[S] = Scattering Matrix 


HN UNES. Twoeport Junction Model of Discontinuity 
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The characteristic impedance of the n 
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mode, Ya’ is der ned 
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where Sen is the kronecker delta. 


BENdefine the coupling coeficients, H 





Whence, 
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+ = ше | 4 3 
On - Al) = = Үй про ДИ 
m m m ES n n n nm 


By representing the A and D amplitudes as vectors, and the 
H coupling coefficients and Y impedance terms as matrices, 


these relationships may be expressed more succinctly as: 


гн JCA” + AT] = [DO + DÍ] 
mn rt n m m 
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J” represents the transpose of rhe matrix [H ]. The 


Tue "LE 
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elements of d may be evaluated by direct substitution of 
the expressions for the normal modes (Eqs. II.A.1 and IT.A.4) 


Мао the definition (II.B.2). From Appendix B we write: 
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Muss scattering matrix for a single junction is defined in 
Appendix C to have matrix elements (which are themselves 


meer ices) as follows: 
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Note fiat eaciseort thesandimidual. elements. H iX T. anda 
themselves matrices. For a single fin-line strip there is 
a scattering matrix associated with the lert boundary 


(minimum 2), Be described above, a transmission matrix 


associated with the width of the strip, [T d, and a 


Scattering matrix associated with the right boundard, 15-1. 


The expression for an equivalent scattering matrix for the 


net effect of LS, 1, [TW], and [Sn], is developed in detail 


in Appendix D. The result is: 
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The equivalent scattering matrix for a fin-line with 
and arbitrary number of fins may be determined by repeate 
Application of equations (11.B.4), with the appropriate 


Fransmission matrices used in the cascade., 


24 





tit. СОШ ТОВ: РТС PROGRAM 


The CAD (computer aided design) program has two functions: 
1) design of a filter (description of physical dimensions) 
based upon user specified characteristics; and 2) analysis 


ОГ а physicaily described fin-line filter., The user has 


Ft 


the option of directly entering the analysis porticn of the 
CAD if filter dimensions have been determined 5v some other 
means. In most cases, however, the user will specify the 
desired response of the fin-line filter (bandedge frequencies, 
passband rippie and upper skirt rejection) and the CAD will 
determine the filter dimensions. These dimensions are then 
automatically passed to the analysis portion of the CAD. 

The design portion of the CAD is built upon the 
theoretical work of Levy [Refs. 6, 7] and contains a modified 
form of a subroutine developed by Shih in [Ref. 2]. The 
design technique is fundamentally a conventional network 
synthesis. The insertion-loss characteristic of the 
distributed filter is approximated by the equation: 


A SA A 
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Е. is the guide wavelength at the filter centeno se abeney 
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>> 
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21052) is the guide wavelength at filter bandedge, ЕСЕ.) 


O 


T. is the first-kind Chebyshev polynomial of degree n 


ть the number © filter Tes Cna tOr: 


The insertion loss equation is evaluated for the input 
specifications to determine the number of resonators required. 
To properly determine the number of resonators for the general 
case, it was necessary to develop subroutines to solve 
simultaneously the transcendental equations for the wave 
numbers of the eigenmodes. The functional dependency of the 
insertion loss equation, above, on guide wavelength is 
EM Ous for the Class II and Class III filter. Thus it 
was also necessary to determine an "effective" guide 
wavelength which accounts for the fact that the wavelength 
changes as the wave propagates from one dielectric region 
into another. It is this effective guide wavelength which 
1s used in Levy's formula. The required reflection 
coefficients of the n+l fins are then derived from the 
Chebyschev polynomial. Next, the scattering matrix for a 
finite septum (Eqs. II.B.3) is evaluated by a bisection 


P C nique to determine the n*l fin dimensions which satisfy 
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these reflection coefficient requirements. Finally, the 
resonator dimensions are determined so that the electrical 
length of each resonator is (1/2) Moog? 
5 

With the dimensions of the filter thus determined, the 
CAD then performs a complete frequency response analysis 
based upon the mode-matching technique described in 
ae ons IIA and II.B. It should be noted that the-equal- 
ripple Chebyshev polynomial expression very accurately 
matches the computed response when the filter is narrow banc 
(less then 4% bandwidth). The agreement between the 
Chebyshev filter approximation and the mode-matching 
scattering analysis deteriorates as the filter bandwidth 
increases. A thorough discussion of this phenomenon is 
given by Shin in [Ref. 8] along with new formulae which 
account for the frequency dependencies of the step-impedances. 
However, excellent agreement between design criteria and 
predicted response has been obtained by employing an 
empirically derived offset to the lower bandedge frequency. 
ME is, if the user specifies the filter Dandedges 1 be 
Бета F2, the design portion is actually conpEEEd 205 


bandedge frequencies: 





The program listing given in Appendix F includes those 
subroutines which were developed as a consequence of the 
previously discussed theory. Subroutine TREQ develops 
eigenmode wave numbers for all three filter classes. 
Subroutine PROPAG determines the corresponding propagation 
numbers. Subroutine SMATRI determines the values of the 
Yo matrix coupling coefficients, and the elements of the 
e tering matrix. Subroutines САЗСАВ, DSUM, and TRANSM 
are used to determine the net effect of several cascaded fins 
and resonators, as discussed in Appendix E. 

Certain subroutines used in the design portion of tne 
program are substantially the work of Shih Rene 
not included in the program listing. Also omitted are 
certain mathematical and algebraic routines required for 


the matrix manipulations. 
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ТУ. PROGRAM VALIDATION 


A. OVERVIEW 

The first phase of program validation was conducted as 
critical subroutines became available. For example, the 
subroutine which determines wave mode numbers was tested 
tor special circumstances for which the values could be 
verified by means of a hand calculator. A series of tests 
were then made to assure convergence to this limiting сасе 
Enoressop Y. C. Shih kindly modified stibroutames s s s 
from his work [Ref. 2] to output values at intermediate 
stages of his analysis routine. These numerical values 
were of interest because they were determined from closed- 
Rom expressions for the reflection сое CE 
employed a residue calculus technique. Thus, the corres- 
ponding values determined by the new program were evaluated 
From a completely different set of expressions and based 
upon a fundamentally different approach to the problem. 
Correspondence between equivalent values was considered 
good to within the accuracy of the numerical process. 

Predicted filter responses were compared to actual 
filter performance to validate further the computer aided 
design and analysis program. Dimensions for Class I and 
Class II filters described by Shih [Ref. 2] and Arndt, 


Radar, |Бет. si were ıhsertsc into ee а па So U DO n CIC 





the CAD and filter performance was predicted over an 


ps 


dpnsopriate fmeguency ramg®e. Tae agreemene between 
predicted response and published E Mer pertormence was 
excellent. Lastly, final validation of the CAD was ассоп- 
plished by means of locally сеет сошенеше сате and 


testing three Class I and two Class il Еше тле filters. 


ERS IL TER CONSTRUCTION AND TEs? 


The locally constructed test filters were designed to 


De resonant within the 8 to 12 GHz band. This frequency 
Ma choice was selected to simpli£y both filter construction 
and test. 

Phe three Class 1 filters were manually CMe with sicissous 
and razor blades from a 2 mil(0.002 ineh) copper foil. Tie 
narrow band Class III filter was made by manually cutting 
EMBy the 2 mil copper layer from 1/8" copper clad substrate 
(Rexolite 1422). The substrate for the wideband Class TIT 
filter was a slab of styrofoam cut from a packing form. 
Conducting fins were cut from 2 mil copper tape and pressed 
onto the styrofoam siabe to form the (ren оделе щите". 

The test fixture used was a specially machined ten inch 
aluminum section of WR-90 compatible waveguide (0.9 by 0.^ 
inch, internal), modeled after a design by Anorr КЕШ 91. 
uide consisted to two longitudinal sections б шше 


This 


Cy 
> 
< 


mere bolted together. The transverse dimension oi the 
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Beer fitter was made equal to Чһе тесе еее ог 
dimension. It was a simple natter, Therefore, o seperete 
dec n ly the test guide halves, slip in the fin-line 
structure, and tighten the clamping bolts. While this 
introduced a 2 mil dimension distortion in the center area 
of the test fixture, no adverse effects were noted. 

The fin-line structure of the Class III filter was 
pressed into the opened waveguide groove. A small spot of 
rubber cement served to hold the dielectric substrate 
against the guide wall. 

Filter electrical performance was measured with the aid 


of an HP-8409 automated network analysis system. 


ни о ЈЕ 15 
A plot of normalized throughput power versus frequency 


lters is shown in Figures IV-1 through IV-7, 


H- 


or seven f 


rh 


ђ 


for Filters #1 through #7, respectively. Filters #1 through 
ің and #7 were fabricated by the authors. For these figures 
mie comtinuous curve is a plot of the Бггасссеа mikt 
response generated by the CAD. The circles contain 
measurement data points. Filters #5 and 76 are representa- 
ihre Filters described by Arndt, et al [Ref. 31. Again, 


the line is the CAD generated plot, but the data points on 


(А) 
Li 
. 


these curves were taken from figures contained in LRef. 


(D 
H 
N 


Tie IV-l contains the physmeal Gesermipuces ot ene Tum 
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Filters #1 and #2 are Class III type structures; filters 
#3, #4 and #7 are Class I; filters #5 and #6 are Class II. 
in general, the shape of the filter response corresponds 
well to the simulated filter. Actual filter passband ripple 
was more than predicted, and varied from a best case of 
ШЕН ПВ (#1) to a worst case of 3.0 dB (#2). Best case 
insertion loss (Filters #1 апа #7) was about 0.3 dB. 
Pegmrications to existing filters and multiple construction 
attempts at the same design suggests that the majority of 
the discrepancy in predicted and actual response was due to 
error in the hand-out filter dimensions. 

The good agreement between predicted and actual filter 
response is readily seen for filters #4, #6 and #7 (Figures 
IV-4, ІУ-6, IV-7). In other instances the agreement is less 
apparent because of the expanded scale of the plot. For 
example, the discrepancy between predicted and actual 
response of filter #1 may be almost entirely accounted for 
|a shift in center frequency of about 150 MHz. This 
150 MHz shirt represents about a 1.7% error in center 
frequency and could result from a construction error in 
resonator dimension of about 0.3 mm. 

Filter #7 is particularly interesting because it 
demonstrates the utility of the program. This fiiter was 
designed to satisfy a requirement of a physics researcher. 


Total time logged on the computer was about 15 minutes. 
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A csv or the output to the Ск ле суеп олет, 

The CAD determined filter dimensions were scribed onto a 
sheet of copper foil and the filter was cut out by hand. 
Total construction time was about one hour. Figure IV-8 

is a copy of the test equipment output, with notations added. 
The minimum measured insertion loss was 0.28 dB at 8.604 GHz, 
Passband ripple was measured at 0.39 dB. The lower 2 dB 
mut was measured at 8.3960 GHZ and the higher 3 dB point 
was measured at 8.823 GHz. Thus, the CAD analysis values 

mem the 3 dB frequency cutoffs of about 8.37 GHz and 

8.85 GHz compare to within 0.3% of the measured values. The 


meeearcher considered the more than 60 dB of out-or-beand 


rejection provided by this filter to be quite adecuate. 


E DISCUSSION 

Experiences related to the design, fabrication and 
test of the fin-line filters motivates us to catalogue a 
number of observations which may be useful to others. 

First, it should be stressed that the filters constructed 
were very crude by industry standards. Dimensional tolerances 
ERE -touHted at * 0.4 mm. thus, the dim@nsional accuracy ОБ 
the smallest fins is roughly bounded by + 40%. The reasonably 
good electrical performance of these filters was, therefore, 
EE the more striking. Also, it should be noted that the 
approximate 150 MHZ shift in center frequency for filter #1 
(narrow band Class III) can be the result of about a 0.2 mm 


Каор in the construction of the resonator dimension. 


SS 








TABLE IV-1 


DESIGN DATA FOR CADOVALIBDENTTON TEST FINEERS 
Filter # l 2 3 i 5 6 7 
a (cm) 1.143 1.143 1.5243 1.143 „7895 21738 1.143 
d(cm) 9.508 02508 0.508 0.508 2912 .00635 0.508 
e, upper) 2a о 1.030 1.000 12008 129 1:59 1.000 
e, (Lower) 12 000 14900 12090 12000 DD 2222 1.000 
Fl (GHz) 8.525 9.950 8450 11,30 15 006 66.0 8.400 
F2 (GHz) 32775 109055 83550 11.05 15960 e 8.800 
 еггірс 3 6 4 4 4 4 2 
* Res 2 5 3 E 3 j 6 
51 (ем) .0959 0.1485 20732 .4995 2.075 . 0347 . 0906 
S2 (cm) .4928 0.7782 „4185. 1.4252 700 . 2088 „4794 
S3 (cm) ¿0959 1.0040 .4185 1.4252 . 7760 . 2038 25961 
S4 (cm) --- 1.0040 ¿0752 „4995 ‚2075 . 0347 .6134 
S5 (cm) --- 0.7782  --- --- --- --- : 5961 
S6 (cm) --- 0.1485 --- --- --- --- .4794 
S7 (cm) --- --- --- --- --- --- 2050905 
Rl (cm) J 5042 Т.1435 1.9991 . 9398 . 8400 .1895 2.0480 
m2 (cm) 1.5044 1.01142 2,0728 ¿9265 .8400 -T1900 2.1227 
RS (беш) mom 1.1104 1.9991 . 93985 .8400 ses 2.1280 
R4 (cm) --- 1.1142 --- --- --- --- 2.1209 
R5 (cm) --- 1.1435 --- --- --- --- 2.1227 
R6 (cm) --- --- --- --- --- --- 2.0480 


m = m 
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Secondly, diffieulties associated with thetfabrication of 


{~ 


the Class III filters suggests the need for an “improved 
construction technique. Any new such technique needs to 
address the two major shortcomings of the method employed. 
These are: 

im Difficulty A Yengitudinally alvenine The wo 
separate fin-line sections to maintain symmetry about the 
x=0 axis. 

2) Persistent probiems associated with inadequate 
ser errical contact between the filter's conducting fins and 
the waveguide's broad walls. 

Failure to correct either of these recurring tendencies 
resulted in severely reduced electrical performance. Thus, 
the original motivation for developing Class III filters 
may have been false. 

Thirdly, dielectric substrates should be tapered. For 
example, the narrow band Class III filter (filter #1) was 
originally constricted on rectangular dielectric slabs. 

The entire insertion loss curve of this filter exhibited 
с 12 dB of ripple. The insertion loss data points 

Pin ед аге for the same filter, Dut wath the addition 

of dielectric ramps on ali four slab ends. Passband ripple 
was thus reduced to about 0.1 dB, and ripple along the 
rejection skirts was no longer noticeable. 


Finally, from an applications view pora, 11 15 





correlation to passband ripple. That is, as the ripple is 
decreased, the fin-line strip dimension decreases. Strip 
width may be increased by either specifying a larger ripple, 
or by forcing additional resonators by requiring sharper 

rejection skirts. Wider strip widths are preferrable simply 


because of achieveable mechanical tolerances, and structural 


integrity. 
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Vee CONCEJOS TONS AND RECOMMENDATIONS 


The numerical analysis technique employed by the authors 
will always produce a result. The question is whetner this 
result has any correlation to the physical world. It was 
found that this numerical analysis yields values essentially 
identical to those obtained by Shih [Ref. 2], who employed a 
completely different approach, and evaluated entirely dif- 
ferent expressions. This suggests that both approaches are 
valid. Furthermore, the good agreement between all 
predicted filter responses and filter performance strongly 
Suggests that the procedure applied was theoretically sound, 
БИШ Әгорегіу executed. 

Mme CAD program produced by this study 15 not am end in 
itself, but a design tool. It should be used to answer a 
number of design questions. It is recommended that an 
organized study be conducted to determine the following: 

1) Under what conditions is the center-line structure 
Beem@errable to a bilateral structure? 

2) For the bilateral fin-line, what determines the 
optimal transverse fin-line location? 

3) Are there optimal values for the relative permittivities 
in the Class II and Class III filter designs? 

4) How sensitive is filter reponse to dimensional 


tolerances? 


u5 





5) Does the current filter design meet engineering 
requirements, or is the use of a design/response optimization 
routine justified? 

5) To what extent can the number of modes retained in 
the numerical analysis be reduced without adversely 
reducing accuracy? 

Finally, to determine more accurately the degree of 
correspondence between simulated and phvsical filter 
performance it will be necessary to construct higher quality 
devices. In particular, the fin-line structure should be 
produced by photo-etching techniques to minimize 


dimensional errors. 


u6 





APPENDIX A 


REGION I EIGENMODES 


m 
— 


The solution to the general wave equation for an EM 


rj 


l 


wave in a guide with disimilar dielectric zones, as shown 


I 


in Fig. II-1 (Region 1), is found in a manner similar to 


that for a simple waveguide. However, in this instance, 


— 


the eigenvalues for the propagation numbers are not 


transcendental 


F1 


y Akbressinte. Instead A Dalir о 


(b 
Уч 
Ж?) 
м 
ШЕ; 
O 
pe 
ct 
к 


equations is found (A.1 and A.2) which must be solved 
Simultaneously to determine each eigenvalue. 


Bor inciaent "EE 4 (m odd) modes, the Z component of the 


1: 


Fi 


шестесіс field (H_) in the wave guide wili have a value o 


сї 
О 


Z 
Zero along the x=0 piane. It is convenient, therefore, 
impose the condition of a "magnetic wall" along x=0, on 


“nich a boundary condition is that 


By referring to the geometry of the waveguide depicted 


| 


ЕЕ. 11-2, we can write the rorm of the E. fields as: 


ик c Em cos ж 
im n J m m im 

Ф = C Sin k "Z TIDER e ik X 
um n um n u 
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At the magnetic wall, x=0, the boundary condition imposed i 
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But since = = 0 for all fields, we can write: 
Ex ула Dn Ecos 415727015075 : 
ке, La m lm dmm Imre 


Whence, A_=0. 
m 


1 


The second boundary condition to be satisfied is that: 


Шс л лла + Тебе колм О 5 
m um m m 


Whence, 
Mm = -C tank а, 
m um 


m 


Therefore, the E fields. :9esxon I is суеп: 
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Or, employing a trigonometric identity, Eyu may be expressed 
as: 
sin k (a-x) 

! = ШІЛ p ° 

ф С -------Ң------- =: R? sink (а-х). 

um m COS k d m um 

um 

Finally, El = Eu at the boundary of the two zones. That is: 


Ecos «and = © sena Fred). 
m ІН! um 


Beem which: 


Likewise, the 4 component of the magnetic field must be 


@eme2nuous across the boundary. So that: 


H4 
ii 
ы; 
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Whence: 


В х ce ug MSS Саа 
m im im m um 


БЕ иле for B above, we find that: 


k Eak d =k cot k (а-а). Сат) 
im im um um 


The second equation which the propagation numbers must 
satisfy is the usual dispersion relation for an EM wave 


propagation through a dielectric discontinuity [Ref. 5]. 


k SENS - Kk СЕ = £). (A.2) 
14 


The amplitudes of the eigenmodes may be determined by 
applying the orthogonality property of the modes. That is, 
trigonometric identities may be used to rewrite the 


expressions for the eigenmodes as follows: 


ien" В. cos Kı „¥ = A, cos Km Due 


ф AS UI XN 
um m um 
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The orthogonalityofthese modes is exploited by integrating 


the product of mode pairs across the width of the guide. 
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E um UNES 7 un 
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APPENDIX B 


ШЕВ ВИШОМ Of TE TEES SIONS TOR HEL TIT ELEMENTS 


The elements of the H-Matrix are defined by: 


а 
Hn я 5 m ex 
where 
ы CEO NES Dc EG 
In — 
= (TI.A.#) 
б = п кш) п а 
un — — 
О КО | 
T OS Km? ОЕ а 
W = { (II.A.5) 
D m 
( ——) OQ NM E) GN NC Nu 
a-d um _ — 
We have therefore: 
d 2 142 
а = ae A, cos ko (cos k, 2) dx 
a 2 ee? 
TO (p sim cM EU E) (sin k (a-x))dx 
а п un a-d um 


32 





142 sin(k,,-k, 29 sin(k, z *k, 4 


2 
Ho: ACD 2  — — Se 
mn d 2 ЕЕ 2 Xy 7 Kin 


2 SUC tky ) (a-d) sintk n Kum’ 270) 





¬ ma o e м в 
un wm 
2 ui? Kn 
H = A (=) (шш COS ап--------1 
mn ud in k 2 aT 
о “Lm 
(шд 372) 
2 1/2 m m 
p cue C=) Ban аа A NA . 
m a-d un 3 2 і 2 
k - Ka 
un Um 





APE Sax aC 


JUNCTION SCATTERING MATRIX 


Consider the waveguide discontinuity as a two-port junc- 
tion problem, as represented schematically in Fig. II-3, 


- + 

where (a ^, a, ) and' (b, 5 ы. ) are the amplitudes of the 
incident and scattered modal fields in Region I and II, 
respectively. 


The electric field is a single-valued function of Z. 


Hence, at the junction boundary, Z=0, we can write: 


Се) 
Early, the relationship for the x-component of the 
magnetic field, Н> tS: 
Ў i = mm E E ba 
n= 1 B E nsl * 
с + со E 
= ш. A ы Do. KO 
же n n n se n n n 
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Toutilize the orthonormality of the eigenmodes, we multiply 


Y 0 | e | 
Et by = and integrate across the waveguide. 
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E ox - т Ж. E 

2 a = а f % pdx = E b + р ЈЕ Шах 
n=l ic О De m Lm ш LEE. E 

И кез з сун ое У (3) 


We define the values H x as: 
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= : ! Д 
Н | иг Y dx 


Likewise, multiplying equation C.2 by m and integrating 


yields: 


The wave impedances of the normal modes, Y may be repre- 
sented as diagonal matrices, the values Hon as an m by n 
matrix, and the amplitudes or the normal modes as vectors. 


Equations C.3 and C.4 may then be expressed as follows 


(subscripts have been omitted for simplicity). 
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EG Е ИО S MN (c.u.1) 


| is solved foe [dp]: 
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ENS tituting [5'] into C.5.1: 
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Whence, 
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ИЕТ shy say ad о IM 
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Ира се, C.4.1 may be solved for La |] and this substituted 
memo C.3.1 to yield: 
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But from standard scattering matrix notation with the 


values denoted as in Fig. II-3, we can also write: 
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IE ODIUM C D 


ЭСАГТЕК ТЕ MATRIX FOR A FINITE CRNO SEPPUN 


A generalized scattering matrix technique in e manner 
БРАШ аг (о that of Shih in Ref. 31 is applied to obtain the 
two port scattering matrix for the finite septum. 

Consider a TE o -type wave incident from Region I at the 
junction A. The field at this junction is comprised of the 
vector sum of all the eigenmodes in each region. Some of 
p -Nenergy will be reflected at the boundary, propagating 
im E) modal fields in the minus Z direction. Phe оюл Ек =з 
Рене епеусу will be transmitted into Region II in e modal 
felds in the positive Z direction. After travelling a 
distance D, the transmitted modes encounter the second 
junction, B, where again some energy is reflected, and some 
transmitted., Multiple reflections of the modes between 
junctions A and B is inherently accounted for in the matrix 
manipulations which yield the scattering matrix for the 
pire septum. 

In C a single junction was analyzed as a two-port scat- 
tering matrix, 15,1. Each of the four elements of this 
Scattering matrix, [S T MS je ES J Jama ES 


Pes A12 Be 


a matrix of infinite dimensions, corresponding to the 


A221: is 


infinite set of eigenmodes in both regions (equation C.8). 
Mie junction B is identical to A, excep. for COTIN IE Ch 


ШО y De shown, therefore, that the scattering matrix for 





the isolated junction B is related to the isolated junction 


at A by the following: 


ES [es 


All po 


[5,12] = [521] · 


The effects of the wave propagating (for propagating modes) 
or attenuating (for evanescent modes) along the strip 
length are accounted for by a transmission matrix, (7,1, 


defined as: 


EL [0 


D 
Co] = 


Fre [I] is the identity matrix, [9] is tae null, 
888 (ТІ іс а diagonal matrix of inrimite size. The 


diagonal elements of LT] are: 
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The combined effect of the single junction A Followed by the 
e A ssion along the strip for a distance D is given by 


the scattering matrix 15,1. 


Vo 1 [S.J CT! 


Pinally, the net effect of the finite septum is determined 


ae follows. 


= 5411 a + 5,12 b Cm) 
Ee = dbi @#2) 
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We define an equivalent matrix for the finite septum, 15.1, 


E 
ban: 
m + + = 
a = 5-11 а + 5-12 а (b o) 
47 = 80) а + 86) 1 (9.8) 
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Wath the equalities given in equations D.5 and D.6, we 


ENDStituce equation D.2 into D.3 to yield: 
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= [Sao JET г 2. 5,111 [3554 d а 
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Likewise, we substitute equation О. into D.i to yield: 
B ! ! ! p г ! 1 
a - 15,4, % 5,4211 - 5і11 3422 35811 Saaı- 
t А t -1 + 
+ 15412111 - $511 S 22 [S344 ay ©: (Bald) 


The matrices which are the elements of the equivalent 
scattering matrix may be written down by simply comparing 


equations D.9 and D.i0 with equations D.7 and D.8. 


! 1 t -1 t 
E С dM | 
321; 5 15,411%15,42 ПІ - 5511 5,221 1534) 5,114 
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3212 5 3445411 - 8444 8455 5312) 
1 ui 1 -1 
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ENDE: É 


CASCADING SEVERAL PWUTE LEGON SEIS 


The fin-line filter is simply a series of finite length 
septa separated by transmission line sections. The net 
effect of the entire filter structure may be represented in 
s ole "toral® scattering matrix which accounts for the 
combined effects of each septa/transmission line section. 
This total scattering matrix is obtained by combining each 
septum equivalent matrix in a manne» precisely analogous 
to that used in D to develop the equivalent matrix for a 
finite septum. That is, in D, expressions for the net 
effect of two discontinuities separated by a transmission 
line section were obtained (equation D.11). It is now 
necessary to express the net effect of two septa separated 
by a transmission line section. Consequently, the equiva- 
lent scattering matrix of two separated septa may be found 
by use of the equation D.11, with the appropriate 
scattering matrices for the septa substituted for the 
Sunetion scattering matrix. More than two septa are 
ar paded by repeated application of D.11. Due care must be 
exercised to assure that the propagation numbers used in 
the transmission matrix correspond tc the appropriate wave 


number and region of waveguide propagation. 
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